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Abstract 

On 

In the multiuser downlink, power allocation for linear precoders that minimize the sum of mean 
squared errors under a sum power constraint is a non-convex problem. Many existing algorithms solve an 
equivalent convex problem in the virtual uplink and apply a transformation based on uplink-downlink 
duality to find a downlink solution. In this letter, we analyze the optimality criteria for the power 
allocation subproblem in the virtual uplink, and demonstrate that the optimal solution leads to identical 
power allocations in the downlink and virtual uplink. We thus extend the known duality results and, 



(N 
> . 

importantly, simplify the existing algorithms used for iterative transceiver design. 
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I. Introduction 

In the multiuser multiple-input, multiple-output (MIMO) downlink, linear transmit/receive 
processing to minimize the sum of mean squared errors (sum-MSE) under a sum power con- 
straint is a well-studied problem. When formulated as a precoder design problem, with implicit 
minimum-MSE (MMSE) receive matrices, the sum-MSE is a non-convex function of the down- 
link precoders. A standard approach used to solve the problem is to find an equivalent formulation 
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in the virtual uplink, wherein the roles of transmitter and receiver are exchanged Q]|-@ . In the 
virtual uplink, the receiver is the Wiener filter and the power allocation subproblem is convex. 

The equivalence of the downlink and virtual uplink problems are enabled by an uplink- 
downlink duality result for the MSE of each data stream. Duality results for linear precoding 
systems were first presented for signal-to-interference-plus-noise ratios (SINR) in @ with single- 
antenna receivers. This work was later extended to MSEs and systems with multiple receive 
antennas in |[T], and subsequently generalized in [6l . The algorithms we focus on are based on 
iterating between the downlink and virtual uplink. Crucially, a common feature is a transformation 
of the resulting power allocation in the virtual uplink to the downlink, while achieving the same 
MSE in each stream in both systems. This transformation requires solving a matrix equation in 
each iteration. 

In this letter, we use the Karush-Kuhn-Tucker (KKT) conditions for the power allocation 
subproblem in the virtual uplink to show that at the optimal point, the powers allocated to 
each data stream in both the downlink and virtual uplink are identical. This result extends the 
known dualities in the multiuser MEMO case. Importantly, this also eliminates the need for the 
uplink-to-downlink power transformation, significantly simplifying the algorithms used. 

This paper is organized as follows: Section HI describes the system model and existing 
algorithms for minimizing the sum-MSE using uplink-downlink duality. In Section [Till we present 
the KKT conditions for the virtual uplink power allocation subproblem, and use the resulting 
expressions to prove the equality of the downlink and virtual uplink power allocations. Section HVl 
wraps up the paper with some conclusions. 

II. Background 

A. System Model with Linear Precoding 

In the linear precoding system, illustrated in Fig. [H a base station with M antennas transmits 
to K decentralized mobile users over flat wireless channels; user k has N k receive antennas. 
The channel between the transmitter and user k is represented by the N k x M matrix H^, 
and the overall N x M channel matrix is H H , with H = [Hi, . . . , H#]. User k receives L k 
data symbols x fc = [xf-i, . . . ,XkL k ] T from the base station, and the vector x = [xf, . . . ,x^-] T 
comprises independent symbols with unit average energy (E [xx ff ] = II (the L x L identity 
matrix), where L = J2k=i ^h)- User fc's data streams are precoded by the M x transmit 
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filter XJ k = [u fcl , . . . , u kLk ], where u k j is the precoding beamformer for stream j of user k with 
Hufcjll = 1. These individual precoders are combined in the M x L global transmitter precoder 
matrix U = [Ui , . . . , IJk] ■ Power is allocated to user fc's data streams in the vector p fc = 
[pfci, . . . ,PkL k ] T an d Pfc = diag [p^]; we define the diagonal downlink power allocation matrix 
as P = diag{ [pf , . . . , p^] }. Based on this model, user k receives a length- iV fc vector yf L = 
HfUv^Px + rifc, where the superscript DL indicates the downlink, and n k ~ £/V(0, cr 2 Ijv fe ) 
consists of zero-mean white Gaussian noise. To estimate its L k symbols x fc , user k applies the 
L k x N k receive filter Vf , yielding the estimated symbols xf L = Vf Hf Ua/Px + Vf n k . 

To minimize the sum-MSE in the multiuser MIMO downlink, we use the virtual uplink, also 
illustrated in Fig. [TJ where each matrix is replaced by its conjugate transpose. In this transformed 
system, we imagine transmissions from mobile user k that propagate via the transpose channel 
H fc to the base station. The transmit and receive filters for user k become V fc and Uf respectively, 
with normalized precoding beamformers; i.e., ||vfcj|| = 1. Power is allocated to user fc's data 
streams as q k = [q kl , . . . , q k L k ] T , with Q k = diag [q fe ] and Q = diag{ [qf, . . . , q£] }. The 
received symbol vector at the base station and the estimated symbol vector for user k are 
Y UL = Eti^iVQl^ + n and x^ L = Eii Uf H^VQ^i + Ufn, respectively, with 
zero-mean white Gaussian noise n ~ CAf(0, u 2 I M ). 



B. Minimum Sum-MSE Multiuser MIMO Linear Precoding 

1 ) Convex Minimum Sum-MSE Precoder Design: The MSE matrix for user k in the downlink 
using arbitrary precoder and decoder matrices can be written as 



H 



(1) 

= Vf J k V k - Vf Hf \J k - Uf H fc V fe + l Lk , 

where 3 k = H^tJPU^H*. + a 2 I, XJ k = Ufc-^/Pfe, and data and noise terms are assumed to 
be independent. The individual MSE terms are minimized using the MMSE receiver V£ = 
Uf HfcJ" 1 . The resulting MMSE matrix is E* DL = I Lk - Uf HfcJ^Hf U fe , and the minimum 
sum-MSE for any choice of U fc is SMSE DL = ^f =1 tr [E£ DL ] . 

The problem of finding the sum-MSE minimizing precoders and power allocations in the 
downlink under a sum power constraint tr [P] < F max is non-convex when MMSE receivers V£ 
are defined as a function of U due to the cross-coupling introduced by the presence of all U, 
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terms in every J k . The authors of (61 demonstrate hidden convexity in the downlink sum-MSE 
minimization problem by optimizing receive matrices (using closed form MMSE precoders) 
and applying a modified cost function. However, this problem is more commonly solved via 
transformation to the virtual uplink, which gives rise to several equivalent problems that can 
be solved using convex optimization. The set of virtual uplink minimum sum-MSE precoders 
and power allocations { (V^, q^J ,k= 1, . . . , K} can be found jointly, by finding the optimum 
co variance matrices = Vf QfcVfc and applying Cholesky or eigen-decomposition [4j. An 
alternative approach finds the optimum precoders V fc and power allocations in an iterative 
manner 10, 01 . The convexity of these problems originates from the decoupling of users in 
the virtual uplink. The MMSE matrix E£ for user k is found using the MMSE receiver 
Uf = Vf Hf J\ 



E 



with Vfc = VkVQk an d J = J2k=i HfcVfcVf Hf + o 2 1m- The resulting minimum sum-MSE is 



*UL 



I Lk - Vf Hf J^HfcV,, 



(2) 
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SMSE 
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tr 



K 



J-^H fe V fc VfHf 



k=l 



(3) 



= L — M + a 2 tr [J' 1 ] , 
which follows from tr [AB] = tr[BA], linearity of the trace operator, and the definition of 
J. Minimizing the sum-MSE thus only requires minimization of tr[J _1 ], which is convex for 
both the power allocation subproblem in q fc and the joint precoder design problem in covariance 
matrices R fc . In this letter, we consider the former optimization problem, which is formally stated 
as 

L 



Ql) 



arg min tr 

91, •••>?£ 



5^h ; hf + cr 2 I M 



. i=i 



(4) 



s.t. qi>0 l = l,...,L, ^^<P max , 



h ll ...,h 1 



Note 



i=i 

where we have defined the effective channel H = [HiVi, . . . , H^V^] 
that the columns in H refer to the effective channel vectors for each individual data stream 
1 = 1 L. 
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2) Uplink-Downlink Duality: The duality results in [0, 01 show that any set of MSEs that are 
achievable in the virtual uplink are also achievable in the downlink under the same sum power 

constraint. This duality can be satisfied by factoring the downlink and virtual uplink receive 

_ i _ i 

beamforming vectors for each data stream I as v z H = p t 2 fyvf and uf = q t 2 Auf ', where 

0i > 0. Any set of feasible MSEs e = diag {ex, . . . , El] can then be achieved by the downlink 

and virtual uplink power allocations 

p = a 2 (e-T>-f3V)- 1 f3 2 l L 

,2/„ T\ fl2,T,\— 1 o2 n 



q=a i ( £ -D-/?^)- i /3^ L , 
where 1 L is the length-L vector consisting of all ones,/? = diag{/3i, 
with 



+ 1, 



(5) 

,/3l}, D is diagonal 
(6) 



where £He [•] denotes the real part of a complex number and 

i = j 



(7) 



III. Equality of Downlink and Uplink Power Allocations 

Based on ©, we see that ^ = \I/ T is a sufficient condition for the equality of p and q. We 
now proceed to prove that this transpose symmetry indeed applies for arbitrary virtual uplink 
precoders V fc as long as the optimum power allocation q* = [g£, . . . , qj] satisfying © and the 
corresponding MMSE receive beamformers are used. 



A. KKT Conditions for MMSE Precoding 

From the objective and constraint functions in ©, the Lagrangian is 



C (q, n) =tr 




(8) 
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and the resulting KKT conditions are 

" hf J 2 h\ 

V£ = - 

h£j 2 h L 



y^g; < fmax, ft > 



(9) 



i=i 

/isum > 0, Hi > 



/'s 



0, fiiqi = 0. 



Here, 0^ is the length-L all-zeroes vector, and is the standard basis vector with a single one 
in the I th position and zeroes elsewhere. The gradient in the stationarity condition follows from 
the identity dJ^ 1 = — J 1 (dJ) J 1 [7 J and the linearity of the trace operator. Thus, 

dtT IJ- 1 } 



dqi 



tr 



dqi 



-tr 



(10) 



-hfj" 2 h, 



B. Conditions for Equality under Optimal Power Allocation 

Having solved © for an arbitrary set of virtual uplink precoders vj, we then find the MMSE 
receive beamformers = J^hiy^qf. With the associated virtual uplink stream MSEs E[ and 
scalars 0i = ■<YqJ\\u*\\, we can then use © to find the downlink power allocation p that achieves 
the same MSEs for each data stream. 

In the case where the optimal power allocation results in one or more inactive streams Si = 
{I e (1, . . . , L) | qf — 0}, this algorithm fails since u* = for / G Sj. However, the same MSEs 
can be achieved for these inactive streams in the downlink by setting pi = 0. The power allocation 
p for the set of active streams £4 = {/ G (1, . . . , L) | qf > 0} can then be found by following 
the specified procedure after deleting the rows and columns from /3, D, and ^ corresponding 
to the inactive streams. 

The coupling matrix ^ is a real matrix whose off-diagonal entries [^L- contain squared 
magnitudes of the end-to-end channel gains from transmitted symbol Xj to the decoded symbol 
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Xj. We observe that \I> = \I> r is satisfied when 



h^u* u* H h 



i j\\ ii i 



or equivalently, 



(12) 



The power allocation terms q* and cancel out, and numerators are equal; thus, an equivalent 
expression for the sufficient condition for p = q is 

hf r 2 \n = hf J-% Vi, j G <S A - (13) 

We rewrite the individual terms in © as hf J~ 2 h; = (/i sum — fj,{). Due to the complementary 
slackness condition (piqi = 0), the dual variables are zero for all active streams I G 1S4 with 
qi > 0. Thus, it follows that 

hf J 2 h, = /i sum VIe5, ; (14) 

that is, (fl"3T) is satisfied, St = and the downlink and virtual uplink power allocations p and 
q that achieve the same minimum sum-MSE are identical. 



C. Discussion 

The equality result presented in Section ITlI-B I was shown to apply for arbitrary vj, as long as 
the optimum power allocation and MMSE receivers are used. It follows that it also applies to the 
optimum covariance-based design, when covariance matrices for each stream are normalized as 
R; = ?*R; and R = Vjvf . This result implies that the virtual uplink to downlink transformation 
stage can be omitted from algorithms using both iterative and joint designs based on a virtual- 
uplink solution Hl-Sl, thus allowing for simplified implementations. 

IV. Conclusions 

In this letter, we have proven that the optimum power allocations for the downlink and virtual 
uplink are identical when minimizing the sum-MSE under a sum power constraint. With this 
proof, we extend the known results in a well studied problem. Importantly, our result simplifies 
existing iterative algorithms, eliminating the solution of a matrix equation in each iteration. 
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